Computational general relativistic force-free electrodynamics

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools one window into the universe's most extreme flows of energy. The discretization Maxwell's equations - needed to make highly magnetized (astro)physical plasma amenable its modeling introduces diffusion. It acts as a source dissipation independent system's physical constituents. Understanding diffusion scientific is key classify their reliability. gives specific limits which results experiments physical. We aim at quantifying characterizing properties our recently developed tool for simulation general relativistic force-free electrodynamics, by calibrating comparing it with other strategies found literature. Our code correctly models smooth waves plasma. evaluate electrodynamics context current sheets tearing mode instabilities. identify that parallel magnetic field ($\mathbf{j}_\parallel$), combination break-down across sheets, impairs resistive find least eight cells per characteristic size interest (e.g. wavelength or transverse width sheet) consistency resistivity origins. High-order allows us provide almost ideal orders convergence (smooth) wave dynamics. layers requires suitable prescriptions sub-grid evolution $\mathbf{j}_\parallel$.